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Abstract 

Machining time is a major performance criterion when it comes to high speed machining. CAM soft- 
ware can help in estimating that time for a given strategy. But in practice, CAM programmed feed rates 
are rarely achieved, especially where complex surface finishing is concerned. This means that machining 
time forecasts are often more than one step removed from reality. The reason behind this is that CAM 
routines do not take either the dynamic performances of the machines or their specific machining toler- 
ances into account. The present article seeks to improve simulation of high speed NC machine dynamic 
behaviour and machining time prediction, offering two models. The first contributes through enhanced 
simulation of 3-axis paths in linear and circular interpolation, taking high speed machine accelerations 
and jerks into account. The second model allows transition passages between blocks to be integrated in 
the simulation by adding in a polynomial transition path that caters for the true machining environment 
tolerances. Models are based on respect for path monitoring. Experimental validation shows the con- 
tribution of polynomial modelling of the transition passage due to the absence of a leap in acceleration. 
Simulation error on the machining time prediction remains below 1%. 
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Nomenclature 



Kinematic and dynamic parameters 





jerk vector 


— V 

A 


acceleration vector 


— V 


feed rate vector 


— V 

X 


position vector 


Jmax,i 


maximum jerk limited by machine 




dynamics on the axis i 


A 


maximum acceleration limited by the 




machine dynamics on the axis i 




acceleration, feed rate and initial po- 




sition on the axis i 


Vf 


programmed feed rate 


VI. 


feed rate reached if Vp is not 




achieved 




feed rate set on the axis i 


Vin 


feed rate of entry into a block 


Vout 


feed rate exiting a block 


Vout 


feed rate exiting a block if Vout is not 




reached 


Vdisc 


maximum feed rate for entering a dis- 




continuity 


Vj 


feed rate entering a discontinuity lim- 




ited by jerk 




feed rate entering a discontinuity lim- 




ited by acceleration 


n 


duration of phase i {ri =Ti — Ti-i) 
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Circular transitions in linear interpolation 



A,0,B 


theoretical programmed path 


A,Q,B 


path described by the machine 


TIT, tol^, 


method to define point Q in accor- 


toly 


dance with the programming method 


R 


radius of the arc inserted on transi- 




tion 


h, h 


length covered before and after the 




transition 




angle formed by segments [AO] and 




[OB] 


dc, jt 


normal acceleration and tangential 




jerk in steady state 


Circular transitions in circular interpolation 


Ri, R2 


radii of circles before and after a cir- 




cle - circle transition 


Parameters used in circular interpolation 





centre of circle 


R 


radius of circle 


Pit) 


current point 


a 


angle covered by the circle arc 


eit) 


current angle 


Ji 


curvilinear jerk limited by the axis i 


Jc 


curvilinear jerk 


Ai 


curvilinear acceleration limited by 




the axis i 


Ac 


curvilinear acceleration 
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Polynomial transitions in linear interpolation 



n frame (o,^,^,~^^ 

A, O, B programmed theoretical path 
{xA^ VAtZa) coordinates of oX in the frame TZ 
{xb, UBtZb) coordinates of Ol^ in the frame TZ 
M point of entry into the discontinuity 

N point of exit from the discontinuity 

Q effective point of passage in the dis- 

continuity 

T time of passage in the discontinuity 

L distance OM 

P{t) current point 

Oi spherical coordinates of the point i 

tolx, toly, tol ztolerance of position on the axes 
Vm , Vn feed rates for entry on M and exit on 
N 



1 Introduction 

High speed machining centres allow for extremely high feed rates to be programmed. However, when 
machining molds or dies, dynamic performances of the machines do not always allow such feed rates to be 
reached. Indeed, according to the quality sought, the segments making up the machining path are often 
extremely short and in such conditions the feed rate reached by the machine will then be limited by the 
NC interpolation time or even the capabilities in jerk or acceleration of the axes [1], [2]. The feed rate 
will then not be constant, leading to considerably lower productivity, a variation in tangential cutting 
forces and impaired quality [3]. Many publications relating to the search to reduce the number of feed 
rate changes base their research on the use of NURBS interpolations. [4] [5] [6] or B-spline [7]. However, 
in the industrial world, linear and circular interpolation remain the most frequently used methods on 
many workpieces. Precise simulation of this type of movement is therefore essential. 

The aim of this work is therefore to propose a comprehensive model intended to simulate the posi- 
tion, feed rate, acceleration and jerk in 3-axis linear and circular interpolation taking the machine/NC 
combination parameters into account. 

At present, NC machine manufacturers [8] propose a displacement law on the axes in trapezoid 
acceleration. This type of command has been studied in the literature by a number of authors writing on 
uniaxial paths with null initial and final feed rates [9] [10]. This movement involves 7 phases (figure 1): 
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Vte [To,Ti] < 



Ai{t) = Ao^i + Jmax,iit — To) 

Vi{t) = Vo,i + Ao,i{t-To) 

~\~2'^max,i ^o) 

Xi{t) = XQ^i + Vo,i{t-To) + 
lAo,i {t - Tof + 

^Jmax^i ^o) 



Vie [Ti,r2] < 



Ji{t) = 

Vi{t) = V^{Ti) + A^a.,,{t-Ti) 
Xi{t) = X,(Ti) + VdTi){t~T-,) 

'^2-^'max,i ^l) 



Jii^^ — J max,i 

A{t) = A,{T2)- Jmax,^{t-T2) 

Vi{t) = Vi{T2) + A{T2)it-T2) 
ytG[T2,T3]{ -\JraaxAi~T2f 

Xi{t) = X,{T2) + V,{T2)it-T2)- 
\A,{T2){t-T2f 

QJmax,i -^2) 



Vie [T3,r4] < 



Vi e [T4,T5] { 



Mt) = 





Ait) = 





Vi{t) = 




, Mt) = 


Xi{Ts) + V,,i (t-Ts) 


Mt) = 





Ai{t) = 
Vi{t) = 

Xi{t) = 



Ai {T4) Jmax,i {t T4) 

Vi (T4) + Ai (T4) [t - T4) 

2^max,i 

Xi (T4) + Vi (T4) (i - T4) 
+lAi (T4) {t - T4)' 



• On phase 6, 



yte[T5,Te] { 



• On phase 7, 



V^{t) 



Mt) 
Mt) 



= 



— A 



- - J rr 



(Tq) + Jmax,i (t — Tq) 
~^'2'^'max.i ~~ 

X, (Te) + V- (Tg) (t - Te) 
+ {Te){t-Tef 



(6) 



(7) 



Figure 1: Control law for a translation axis 
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Acceleration starts by increasing up to its maximum value Amax,i- The slope for this acceleration 
is called jerk. This is noted Jmax,i (phase 1). An acceleration level is then reached with value A,nax.i 
(phase 2). This then diminishes down to zero (phase 3). During the fourth stage, acceleration is null 
and the feed rate attained is equal to the programmed feed rate. Phases 5, 6 and 7 are symmetrical with 
phases 3, 2 and 1. Time resolution of this system of equations allows the times for each of the phases 
to be calculated readily (section 2.1). Parameters Amax,i and Jmax,i are specific to the dynamics of the 
machine used. 

In practice, the blocks follow on from each other and the feed rates at the start and end of the blocks 
are rarely null. Furthermore, according to the length of the segment and the feed rate to be reached, 
maximum acceleration or the programmed feed rate are not always reached. Thus, segments are often 
crossed in which some of the seven phases do not exist. This general case is rarely addressed [9] [11]. 
Modelling of the machine's behaviour taking these aspects into account is proposed in section 2. 

In the case of 3-axis machining, the path of a segment is covered simultaneously on 3 axes. Take, 
for example, the case of a block with displacement on ~t and 7". To ensure the quality of the machined 
workpiece, the path actually followed needs to be monitored in relation to the programmed path. To do 
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so, synchronization of the axes is needed. Indeed, both axes must reach the final position at the same 
instant. Thus, on 3 axes, a feed rate, an acceleration and a maximum jerk have to be calculated for each 
of the axes as a function of the slowest axis [12]. Using the model in 7 phases, the same time for each of 
the phases will be obtained on each of the axes, allowing the path to be followed. 

Circular interpolation is also widely used. However, as far as can be ascertained, there are no publica- 
tions covering changes in feed rate on a circle while also ensuring monitoring of the path to be followed. 
Modelling of such a case will be presented in section 2. 

Moreover, a machining operation can be broken down into a multitude of linear or circular blocks. 
This thus poses the problem of the crossing of transitions between discontinuous blocks tangent to each 
other. From a purely kinematic point of view, the exact passage by the programmed points requires 
precise arrest of the machine at the end of each block. This behaviour is not permitted in practice as 
it implies repeatedly slowing down and thus a loss in productivity. In addition, the jerks are prejudicial 
to the quality of the workpieces manufactured as well as the lifetime of the cutting tools used. Thus, in 
numerical commands, there is a tolerance on pursuit of the path, figure 2. For a programmed path A - 
O - B, the machine will cover the path A ~ Q ~ B. This point Q can be defined in a number of ways: 
by the distance TIT, tolerances tolx and toly on the axes and the distances di or d2- 



Figure 2: Model of circular pass for a discontinuity at a tangent 




In the literature, a circle arc is used to model this transition [14]. Radius R of this circle arc is 
calculated in relation to the tolerance permitted by the NC and lengths li and I2 of the blocks before and 
after the discontinuity: 



cos 



R = min TIT- 



TIT 



1 — cos I 



(8) 



,s(f)'sin(f) 
where I — min (li, I2). 

From a kinematic point of view, the course of a circle at constant feed rate V runs at constant 
centripetal acceleration Oc = Deriving this acceleration, a constant tangential jerk jt = ^ is also 
obtained. As this acceleration Cc and jerk jt are limited on the machine by maximum values AmaxA and 
Jmax,i, for & given radius, the feed rate in the discontinuity will thus be limited: 



Vdisc — niin ^Vp , Va — \/ A„iax,iR, Vj — ^ Jmax,iR^ 



(9) 
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This modelling allows the maximum speed of passage to be expressed as a discontinuity. However, it 
pre-supposes a jump in acceleration: at constant feed rate on segment [AO] acceleration is null (phase 4) 
while at constant speed on a circle, the projection of the acceleration vector on the axes can reach Vp/R. 
This leap in acceleration on crossing the transition is not observed in practice. Thus, a circular model 
cannot be used to simulate precisely the position, feed rate, acceleration and jerk along the discontinuity, 
even if it gives a good approximation of the drop in feed rate. In part 3, a polynomial form of modelling 
for transitions between non-tangent segments will be presented. 

In circular interpolation, the transitions between the circles of different radii need to be modelled. In 
this case, the same type of limitation arises: two tangent circles with different radii are discontinuous in 
curvature. On crossing the discontinuity, there will thus be a jump in acceleration. Pateloup [15] proposes 
a model to determine the minimum feed rate needed to cross the discontinuity taking into account the 
jerk j, the radii i?i and i?2 of the two circular portions, and the interpolation time 5t of the machine: 



This model is also used to cross a transition between a segment of a straight line and a circle arc when 
they are tangent. The model is again taken up in the algorithm proposed by Tapie [1] to calculate the 
entry rate into this type of discontinuity. This method for crossing discontinuities will be validated in 
section 4. 

With the aim of simulating a complete path (path including blocks and transitions), Lavernhe [12] 
proposes a method allowing the machine's dynamic behaviour to be computed using a formalism with 
inverse time. Integration of the NC cycle time in this method allows it to predict the control jerk value 
to be predicted for each of the periods. From this are deduced the plots for acceleration and feed rate. 
Furthermore, his model takes the predictive functions available on NCs into account. The model for 
crossing of discontinuities in tangency is that described previously (circle arc). 

Another solution is to identify the servo-system model for the machine/NC combination [13]. This 
approach appears difficult to implement given the lack of data provided by NC manufacturers. Indeed, 
to apply this approach would require precise knowledge of the slaving flow diagrams for the axes and 
especially the various correctors used. Where appropriate, tests need to be conducted to identify the 
transfer function parameters. 

Finally, a third method involves modelling directly the laws described in the previous sections. The 
difficulty in implementing these models lies in calculating for each block the time for each acceleration 
phase as well as the jerk on each axis. Integration of anticipation is no easy matter. 

To sum up, the paths of linear blocks are clearly described in the literature. However, the general case 
(path of a segment at non-null initial and final feed rates) is not studied. Furthermore, no information is to 
be found on the path of circular blocks. With respect to discontinuities in tangency, the model for passage 
in a circle arc allows the feed rate on crossing the discontinuity to be quantified but does not enable laws 
for feed rates and accelerations to be modelled. The idea is to propose algorithms that, on 3 axes, show 
how to go from a feed rate on block entry V/„ to a feed rate on block exit Vout while attempting to reach 
the programmed rate Vp both in linear interpolation (section 2.1) and circular interpolation (section 2.3). 




(10) 
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A model is then proposed for passage into discontinuities in tangency between two straight hnes (section 
3). Finahy, tests vahdating the simulator are presented. 



2 Modelling NC behaviour in linear and circular interpolation 

In this section, the laws for displacements, feed rates, accelerations and jerks to cover a uniaxial 
segment from a feed rate V/„ to a feed rate Vout passing through a feed rate Vp are modelled. Passage 
in the case of a 3-axis segment is then studied. The results are then adapted to circular interpolation. 



2.1 Modelling uniaxial linear segments in the general case 

In the general case of a toolpath for a segment in 3-axis machining, the feed rates at entry, middle 
and exit of a segment will not be identical. In addition, according to the jumps in feed rate to be crossed 
and the length of the displacement to be made, the feed rate will not necessarily be reached. As a result, 
resolution of equations 1 to 7 allows all existing cases to be identified (figure 3 and appendices). Thus, for 
L, Vp, Vin and Vout given, it can be determined whether Vout and/or Vp will be reached going through 
Amax,i or not. This means the duration for each of the phases can be known. The contour for jerk, 
acceleration, feed rate and position on one axis can thus be retraced. 

Figure 3: Set of cases of resolution 
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Is ArnaT.i rcachablc to reach Vq^^^'' 


YES 


NO 



\Vln.Vout. Vf~ 



Ls Vuu, reachable ? 
Eq. 59 or Eq. 62 or Eq. 03 



Calculation of 1 '' 



Is AmaxA reacha.blo 
to reach VLl 



Is Vf reachable ? 
Eq. 70 or Eq. 72 
or Eq. 74 or Eq. 7ti 



Is Ijuiij.i reachable 
to reach Vpl 









Is A,rn,j:,i reachable 




Im A,„a.i_i reachable 




to reach Vqhj? 




to reach Vovt' 




to reach Vout- 




IV'O.. - Vp\ > ^=ji 




\Vo.i-VFl > ^=ff^ 




\Vo.t-VF\>^^ 


YES 


1 NO 


YES 1 NO 


YES 


1 NO 



la A,„a-i.i readjable 
to reach Voui.'^ 
IVb.i- VfI > 



NO 



Case. 1 I I Casf2 | | Cast- 3 | | Ca« 4 | | Case 3 | | Case C | | Case 7 | | Case 8 | | Case 9 | | Case 10 



2.2 Passage to 3 axes 



The algorithm given in the previous section is valid on a uniaxial displacement. In this case, values 
Amax,i and Jmax,i wiU derive from the machine characteristics. On 3 axes, synchronization will be 
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needed. This means that the times for each of the 7 phases are identical on ah axes. As a resuh, for a 
given displacement at a given programmed feed rate, the NC will recalculate a set feed rate, acceleration 
and jerk for each of the axes. 

The modelling method followed is thus as follows: for a given segment, the displacement to be made 
on each of the axes is calculated. Using the results of the previous section with maximum acceleration 
and maximum jerk, it will thus be possible to determine which of the 3 axes will be the slowest. Then, 
using the results Lavernhe offers [12], feed rates Vi, accelerations Ai and jerks Ji on the axes limited i 
can thus be calculated as a function of the distance Lum to be covered on the limiting axis and distances 
Li to be covered on the limited axes: 

T Li ^ Li , , 

^ — ~p Ai — A^fiQ^^i Ji — — Jrrtax^i v-^-^) 

All the elements used to simulate tool paths in 3-axis linear interpolation have thus been presented. In 
what follows, the case of circular interpolation will be studied. 



2.3 Modelling tool paths defined by circle arcs 



As has been seen, few data are given as to simulation of displacements on a circle. The problem 
involves understanding how the axes of the NC behave to follow a circular path and especially how 
decelerations and accelerations are made when following a circle. To this purpose, the laws governing a 
circular movement can be stated. 

Let P be a point moving over a circle arc of angle 0, centre O and radius R in the frame (e^, e^, et) 
(figure 4). 



Figure 4: Parametrization for a circular path 






B £e 




/ \i P 













A 



This gives: 

0P = Ret (12) 
Note 6 the angular law such that Vt G [0,T], 9{Q) = and 0{T) = a. The feed rate, acceleration and 
jerk vectors can thus be expressed as follows: 



dt 

d^oP 

dt^ 




(13) 
(14) 

(15) 
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When the set feed rate is reached, d6(t)/dt = Vp/R- Nevertheless, the angle law 9{t) is unknown on 
passage from null feed rate to the programmed feed rate. As a result, using a law in 7 phases is proposed, 
like that considered in section 2.1. In these conditions the jerk effect no longer corresponds to a jerk by 
axis, but a "curvilinear" jerk taking into account the influence of 2 or 3 axes along the plane in which 
the circle is made. This jerk is thus the result of several contributing axes. This means the tool path has 
to be expressed as a projection on the machine's axes of translation. This is done through two moves in 
frames. 

Let Ml and M.2 be the respective matrices for passage of the frame i^) towards frame 

{1^ ,lf ,~^) and the frame (ep, e0,et) towards frame (it ,lt ,1^). 

Ai is noted as the matrix for passage of the frame (~af,~^,^) to frame (cr, e^, et): 

M = {MiM2)~^ (16) 



2.4 Calculation of jerk and curvilinear acceleration 

Curvilinear jerk and curvilinear acceleration do not form part of the machine parameters, but they 
are the result of the contribution made by accelerations and jerks for each axis in movement. According 
to the zone in which the circle is completed, one or other of the axes will be limiting. 

These parameters are calculated at the start and the end of movement as that is where they are most 
significant. Indeed, for the path of a circle arc, one needs to switch from a null normal acceleration to 
a normal acceleration equivalent to Vp/R at the start and at the end of the path, which would require 
an infinite jerk at the start or the end of the path. As this can only be considered on certain types of 
machining centres, the maximum jerk that can be reached at the start and end of the path in consideration 
of the characteristics of the axes is thus calculated. 

Consider a circle arc made in the frame {it ,1^ ,1^) from an angle 0(0) = to an angle 9{T) = a, T 
being the overall duration of the path. In this general case, the matrix is a matrix for rotation of 
the orthonormed frame (cr, eg, et) to the orthonormed frame ('5', l/, ~^). This matrix is thus orthogonal 
and can be inverted. Note {ux,Uy,Uz), {vx,Vy,Vz), {wx,Wy,Wz) the coordinates of the vectors it , if , 1^ 
in the frame (lzf,~j/,^). 

Using a model in 7 phases for angular acceleration and assuming initial acceleration to be null, the 
following is obtained on the first phase, Vt G [0,Ti]: 

^(t) = Jet 
' ^ ■ 1 2 

^ e{t)^eint + \Jct^ 

with Jc the curvilinear jerk. Carrying over equation 17 into equation 15 on i = and projecting onto 
the machine axes, the following is obtained: 
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R J. 



4 



^ r(Jc- 01,] ^ 



(18) 



(1^,1/,^) 

Jerk is limited on each of the axes and curvihnear jerk wih depend on 2 or even 3 axes. The latter 
will thus depend on the axis that will be limiting. The following will therefore obtain: 



Ji 



J, 



max,x A3 



RVx 



J2 



J n 



RVy 



J n 



(19) 



Similarly, J4 , J5 , Jg are calculated at the end of movement to obtain: 



J, = min (J,) (20) 

»e[i,6] 

Curvilinear acceleration can now be calculated. Noting 61, the position reached at the end of phase 
1 and 6*1 the feed rate reached at the end of phase 1, the following will obtain in phase 2, V< G [ri,T2]: 



0{t) = 6*1 + Oit + \A^e 
Carrying over equation 21 into equation 14, this gives t ^ Ti: 



0W = o 



df^ 



RAe'4 - ROl^'^r 

— {ux COS 9i + Vx sin 9i ) R9l 
+ {uy cos 9i + Vy sin 61 ) RAc^ ■ li^ 
+ ^ — {~Ux sin 01 + Vx cos6'i) Rdf 
+ {—Uy sin 61 + Vy cos 9i) RAc^ ■ IJ' 
- WxR9l 



WyRA, 



Thus: 



(21) 



(22) 



Ai 
A2 
A3 



COS 9+Vx sin 6)R6^ 
{Uy cos 6-\-Vy sin 6)R 
r,y + ( — sin Bi+V:c COS ei)R.e\ 
{ — Uy sin 6i-\-Vy cos 6i)R 



WyR 



The same method is adopted for the deceleration phase, giving: 



(23) 



A, = min {A,) (24) 

To conclude, in the modelling proposed in circular interpolation, the angular law follows a movement 
in seven phases, with maximum jerk being the calculated value Jc and maximum acceleration being the 
calculated value Ac- Knowing all the parameters of a circle arc (point of departure, point of arrival, 
radius, etc.), the machine's dynamic behaviour can be simulated throughout the circular path. 
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This part allows the laws of position, feed rate, acceleration and jerk on unique blocks in circular 
and linear interpolation to be simulated. What remains is to model the junctions between blocks. It has 
already been shown that the model for transition between two tangent paths functions. The following 
part of the article will cover how to model the passage between two linear segments. 

3 Modelling of the transition between two rectilinear blocks 

In the literature, transitions between two segments are always modelled by circle arcs, though this 
does not seem to match the real behaviour of the machine. Knowing that NCs are capable of describing 
polynomials of degree 5, it is suggested that they be used to model discontinuities in tangency. This 
model should allow a criterion for feed rate for entry into the discontinuity to be defined as also the 
contour for position, feed rate, acceleration and jerk in the discontinuity. 

3.1 Notations and hypotheses 

Consider the paths of two consecutive non-aligned segments AO and OB (figure 5). 

Figure 5: Parametrization of passage of discontinuity 

z 










^ '"^^^^N 





The points M and iV can be expressed as follows in this frame: 



OM = 



sm cos Om 



\ 



~ots 



L sin LpM sin 0^ 

V L cos ipM J ^ 

( L sin ip NCOS 9 N \ 
L sin (fN sin^AT 

V L cos ipN J ^ 



= Lit 



Lit 



(25) 



These notations allow the working assumptions to be determined. Firstly, it shall be considered that: 

IIOMII ^ ||OiV|| L (26) 
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Then, considering the problem to be symmetrical, this gives: 



OQ ^ \\Oi 



IZ" + 17" 



The coordinates of point Q can then be expressed in the frame TZ 



Qy 



(27) 



0(^ - \\0(^\ 



( sin 



cosf 



sm LpQ sm f 

C0S(y5Q 

arccos ■ 



V 



n 



Si Wx > , 0Q ~ arcsin ■ 



(28) 



Si Wx < , 9q = TT — arcsin 



The direction of the vector 0(^ is thus fully determined by the vectors it and it. Its norm now needs 
to be determined. This is done by the tolerance granted the machine on passage of the discontinuity 
(section 1 and figure 1). In this instance, tolx, toly and tolz denote the maximum tolerances for passage 
used by the NC on each of the axes and ^, which means that: 



-JI < tolx 
0(^-^ < toly 
< toh 

The norm of vector 0<^ can thus be calculated as follows: 



(29) 



IOC 



toK 



toL, 



toL 



(30) 



I sm ipQ cos tiQ I I sm <^q sm t/g | | cos ^pg I 
Finally, due to the problem's symmetry, it can be considered that the entry feed rate in the disconti- 
nuity ||Vm|| and the exit feed rate ||Viv|| will be equal: 



Vm 



(31) 



As points A, M and O, as well as points B, N and O are aligned, angles lpm, 9m, 'Pn and 0^ can 
thus be calculated: 



ZA 



LpM — arccos ■ 
Si > , 9m = arcsin , „ 
Si < , 9m = tt — arcsin , ^ 

Si xb>^ , = arcsin , , 
Sia;B<0,(?Ar=7r — arcsin . 



(32) 
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3.2 Formation of the equation 

Using a polynomial representation, the equation for the position, feed rate, acceleration and jerk of 
the point P in the discontinuity thus takes the following form, Vi G [0,T]: 



doP{t) 
Jt 



df^ 



' v^5 




y' 






-Q Cit 


. . c 

Et 




Et 








2-^1=2 




Ei=2 


l<. U fi 

(■i-2)!"*'' 


2^1=2 




Ei=3 


(■i-3)!"*'' 




(i-3)!"*'' 


Z^i=3 


(i-3)!^*'' 



i-2 
i-2 
i-2 

i-3 



(33) 



(34) 



(35) 



(36) 



3.3 Boundary conditions 

In order to resolve this system, eight boundary conditions are used: 




Equations 37 and 38 translate entry into the discontinuity. Equation 39 translates the problem's 
symmetry, meaning that the point parametrized by tolerances of passage is reached half way through the 
time of the path. Considering that acceleration is null at entry and exit of the discontinuity, equations 
42 and 43 are obtained. To conclude, equation 44 translates symmetry of the acceleration contour. 
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3.4 Resolution 

This thus involves resolving a system of 24 scalar equations (projection of equations 37 onto 44) whose 
unknowns are: 

• 18 coefhcients Oj, 5,;, Ci of polynomials, 

• the norm L of vectors OA^ and OI^, 

• time T for passage of the transition. 
Resolving the system gives: 



16Qx sin cos 9m 

3 (sin ipiv cos ^^v+sin ip m cos 9m) 

16Qy sin sin ^Aj- 

3(sin c/jjY sin 0jv -f sin ifiM sin ^^a7 ) 

16Qz cosipAf 

3(cos ipN+cos ipAf ) 

ai = - Vin sin ipM cos 6'm 
6i = - Vin sin (^M sin 0m 

Ci = -Vin cos (fiM 



(45) 



(46) 



(47) 




T 



( l\ 

L 

V i / 



9 (sin cpiv cos ^iv+sin ipA/ cos 9m)^VP^ 
1024Q2 

9 (sin c^jv sin +sin tpA'/ sin ^a/)^ V^^^i 
1024QJ 

9 (cos +COS Pm)^ 

1024Q2 

27(sin i^iv cos S^v+sin a/ cos 9 m)^Vi 
65536Q'' 

27(sin 

65536Q:^ 

27(cos if iv +COS if M )*^ V'/t 



sin cpiv sin f?jv+sin (^a/ sin 9m)'^V^^ 



""65536QI" 




OP{T) = 
/ 32^ 



3V7„ (sin ip^ cos ^^^v+sin ipM cos ^a/) 
32Q„ 



3V7Ti(sin (piv sin ^iv+sin c/^m sin 
\ 32Q^ I 

\ 3V/Ti(cOS93iV+COSV?M) / 

Oi^(Z)^0^ ^ 

16Qx 

3(sin c^jv COS ^jv+sin ^m cos (?a/) 

16Qy 

3 (sin ipj<i sin f?jv+sin ^m sin t^Af ) 
3(cos:^iv+cos(^A-/) 



(48) 



(49) 



(50) 



(51) 



(52) 
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Due to the relation between Q^, Qy and Qz, the 3 expressions allowing L or T to be computed give 
the same results. With respect to the feed rate on entry of the discontinuity V7„, it is ex ante equal to 
the feed rate at the end of the upstream block. Nevertheless, it can perhaps be limited by maximum jerk, 
maximum acceleration and the length of the blocks upstream and downstream from the transition. This 
feed rate can now be calculated. 



3.5 Feed rate on entry in the discontinuity 

Expressing the problem as an equation allows the entry feed rate to be calculated when it is limited 
by acceleration or maximum jerk: the assumption is made that maximum jerk is at point M and that 
acceleration will be at its maximum at Q. Care must therefore be taken to ensure that the machine's 
capabilities are not exceeded at these points. 

If the maximum jerk is reached on each of the axes, the equation gives: 



^(0) = J„ 



a,] 




Jmax.x \ 






J max.y 


V C3 ) 




^ Jmax.z J 



(53) 



According to the case, resolution of this system allows the entry rate limited by an axial jerk to be 
determined using the results given by equation 48: 

^lim,j 3|sin c^jv cos f^jvH-sin ipA/ cos 1 min ^ (54) 



Similarly, it is known that on the discontinuity, acceleration is at its maximum on T/2. The entry 
feed rates that will be limited by acceleration can then be determined by resolving the following equation: 

By inserting this condition into equation 35 and using the results given by equations 48 and 49, the 
following is obtained: 

Viim a — "^Tv- o — 1 — n — r min ( 

3|sm i/Piv cos tfjv+sm i/^M cos t^A/| \^ (56) 

y/ Q x-^max ,x ; \/ Qy-^rnax,y ; ^^ Q z-^max,z^ 

Thus, maximum acceleration on each of the axes also leads to limitations on the feed rate on entering 
the block: 

Vin = min (Vp, VHjn,j,Vlim,a) (57) 

Finally, a last case remains: the transition length calculated L can be greater than the length of 
the segment upstream or downstream. An additional condition is thus imposed: if L > ||074||/2 or 
L > ||Oi^||/2, the coordinates of point Q are recalculated taking L ~ min |^||Oy4||/2, ||Oi^||/2^ from 
equation 52. Equations 54 and 56 then allow the limit feed rate on entry in the discontinuity to be 
recalculated. 

The proposed model is now complete. Experimental validation is proposed in the following section. 
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4 Experimental validation 

Measurements of position profiles, feed rates and acceleration were made on a machine with the aim 
of validating modelling. Tests were conducted on a DMU 50 eVo 5-axis machining centre equipped with 
a Siemens 840D Numerical Control. The characteristics of the NC and machine combination are given in 
table 1. This NC allows measurements of position, feed rate and acceleration to be made for a maximum 
period of 10 seconds. 

Table 1: Characteristics DMU 50 eVo - Siemens 840 D 



Dynamic characteristics of axes 



Maximum feed rate 




- 50 m/min 


Maximum acceleration 


A 


= 9,g 


1 m/s^ 


Maximum jerk in translation 


J max^i 


= 40 


m/ 


Maximum jerk on passage of 


J ^raax 


= 60 


m/s^ 


a discontinuity in curvature 








Characteristics of the NC 


Tolerances on the axes 


Tol = 


= 0,01 


mm 


Interpolation time of NC 




2ms 





A test campaign enabled the proposed models to be validated. In particular, feed rates ranging from 
500 to 10000 mm/min and machining tolerances from 0.1 to 0.01 mm on various paths were tested. More- 
over, as the NC allowed certain dynamic parameters to be modified, jerk was also subjected to variation. 
In the present publication, only two tests are presented (figure 6^): one test in linear interpolation and a 
second one in circular interpolation. The first test sought to validate the simulation method adopted as 
well as the method for passage of discontinuities using polynomials. The second confirmed the modelling 
of behavioural laws on circular portions. The programmed feed rate was 5000 mm/min on the cases 
studied. 

The cutter - workpiece feed rate profiles measured and simulated were plotted on figure 7. What is 
immediately striking is the model's faithfulness to the experimental curves. 

In linear interpolation, the path time calculated is 4237 ms and the time measured on the machine is 
4232 ms i.e. an error of less than 0.2%. Error between simulation and measurement reaches a maximum 
level on drops in feed rate at passages between discontinuities (table 2). Although errors between the 
rate measured and the rate simulated in the discontinuity can appear to be significant, such error on the 
drop can be seen to represent less than 1%, the profiles are seen to be closed and error on the time of 
pass through a discontinuity is low. 

Figure 8 shows the plots of positions, feed rates and accelerations in projection on the axis. An 
excellent match can be seen between the model and the measurement. Figure 9 represents a zoom onto 
the zone on crossing the discontinuity between the 600th and 900th milliseconds to pass points 2-3-4. 
Note that on the measurements, the feed rate diminishes on the segment to reach a rate of 281 mm/min. 

^In the case of the circular interpolation, X and Y are the coordinates of the final point 
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Figure 6: Simulated paths 



Z (mm) 



Y (mm) 




X(mm) 



X 


Y 


z 




50 





100 




40 


10 


95 




60 


10 


90 


mm) 


70 


-20 


85 


> ' ■ 


40 


-30 


80 




20 


20 


85 




60 


40 


90 




90 





95 




110 


30 


100 






X(mm) 



X 


Y 


R 


-5.0 


10.0 


7.5 


5.0 


30.0 


25 


30.0 


15.0 


15.74 


15.0 


-5.0 


26.56 


-10.0 


-5.0 


29.51 


-20.0 


5.0 


20.73 


40.0 


20.0 


31.50 


35.0 


5.0 


21.62 



3- axis paths in hnear interpolation 



2-axis paths in circular interpolation 



Figure 7: Measured and simulated contours for cutter - workpiece feed rates 
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The calculated feed rate is 263 mm/min, representing an error of 6.8%. The discontinuity is then entered. 
The feed rate measured at the end of the discontinuity is 91 mm/min. The rate calculated is 85 mm/min 
i.e. 7% of error. Finally, observation of measurement of the acceleration curve shows that there is no leap 
in acceleration on entering the discontinuity as would tend to suggest passage of a transition by a circle 
arc. This corroborates the validity of using the model in which acceleration remains continuous. 

In circular interpolation, the radius values chosen only highlight a drop in feed rate on passage of the 
first discontinuity. Error between the calculated value and the simulated value is here only 2%. Error on 
the difference between the total simulated time and the total calculated time is here also less than 1%. 
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Table 2: Errors calculated on the cutter - workpiece feed rate on passages of discontinuities (Feed rate in 
mm/niin) 

Feed rate in the transition Drop in feed rate 





Measurement 


Simulation 


Error 


Measurement 


Simulation 


Error 


1 


71 


72 


1,4 % 


4929 


4928 


0,02 % 


2 


256 


222 


13,28 % 


4744 


4778 


0,7% 


3 


150 


135 


10,67 % 


4850 


4865 


0,31 % 


4 


193 


148 


23,31 % 


4807 


4852 


0,94 % 


5 


161 


150 


6,83 % 


4839 


4850 


0,23 % 


6 


196 


158 


19,39 % 


4804 


4842 


0,79 % 


7 


126 


91 


27,78 % 


4874 


4909 


0,72 % 



Figure 10 shows positions, feed rates and accelerations measured and simulated on the spiral pro- 
grammed in circular interpolation, in projection on the ~i axis. The acceleration profile perfectly illus- 
trates the leaps on passage of discontinuities in curvature. The Pateloup model used for simulation is 
thus fully confirmed. 

Figure 11 represents a zoom on a passage through discontinuity. The model can be seen to faithfully 
reflect the measurement. 

5 Conclusion 

The present publication proposes a model allowing jumps of discontinuity in tangency to be overcome 
together with a model to simulate the axial behaviour of a machine in circular interpolation. This 
means that the position, feed rate, acceleration and jerk on each of the axes catering for the dynamic 
parameters of the machine/NC combinations can be simulated on 3 axes and cutter - workpiece contours 
be determined. This gives faithful simulation of the machine's dynamic behaviour in 3-axis machining in 
the commonest programming modes. 

Moreover, the present study stresses the problem of using linear interpolation when finishing complex 
workpieces. In this context, programmed tolerances are often extremely low and segments thus extremely 
short, leading to repeated reductions in feed rate. This problem is only very partially resolved by NCs 
using predictive functions. 

This means that beyond a certain segment size, the type of programming needs to be changed. Future 
works will therefore study the machine dynamics in polynomial interpolation and integration of this type 
of programming into the simulator. The possibility of integrating the behaviour of series 5-axis machines 
into the model can then be considered. Finally, in the longer term, the simulator should allow for 
machining paths to be optimized so as to boost the productivity of ultra-high speed machines. 
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Figure 8; Comparison of contours measured and simulated on a rectilinear path 



no 

100 
90 
80 
70 
GO 
50 
40 
30 
20 











Siimilatio 


1 

e 



















































































































































500 1000 1500 2000 2500 3000 3500 4000 
Time (ms) 

Position 




500 1000 1500 2000 2500 3000 3500 4000 
Time (ttts) 

Feed rate 



2000 
1500 
1000 
500 


-500 
-1000 
-1500 
-2000 











Simulation 












































f 




















T 










1 



































500 1000 1500 2000 2500 3000 3500 4000 
Time (ms) 

Acceleration 



Acknowledgements 

This work was carried out within the context of the working group Manufacturing 21 which brings 
together 11 French research laboratories. The topics addressed are as follows: 

• modelling of the manufacturing process, 

• virtual machining, 



emerging manufacturing methods. 



21 



Figure 9: Passage of a discontinuity at a tangent 
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Figure 10: Comparison of contours measured and simulated on the circular path 
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Figure 11: Passage of a discontinuity in a curvature 
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Appendices: General case for the path of a uniaxial segment 

The non-linear system described in section 1 therefore needs to be resolved. This resolution involves 
resolving the 10 cases shown in figure 3. 

The assumption is made that acceleration is null on entry and exit of the block. This is verified when 
the segments are long enough to reach the feed rate programmed at the end of the segment (section 3.5). 

Resolution of this system means the times for each of the phases can be determined so as to plot the 
cutter - workpiece feed rate profile. In some cases, the programmed feed rate will not be reached and the 
rate actually reached will be sought. 

According to the distance to be covered and the jumps in feed rate to be crossed, 10 cases were listed. 
From resolution of equations 1 to 7, the following obtains: 



Case 1: Vout cannot be reached. Vp^f is reached by reaching Amax^i (figure 12): 



< 



\yout-yin\ > 

2Vlr 



^ < L < 



J™ax,.(v^„t-y/j+A^„,_,(y,„+vb„t) 



2A„ 



iJrr 



(58) 



V 



• (yAraax,iJmax,i^ + ^^max. 
~'^-^'max.,i'^rnax,iVln + 



■V, 



1/2 



(59) 



n = T3 = 



T2 



T4: ^ n = Te = Tj = 

• Case 2: Vout cannot be reached. Fq^j is reached without reaching A. 



(60) 



(figure 12): 



\Vout-Vln\ > 



+ 



(61) 



\Vout-Vi^\ < jff fyf- 
L < {Vout + Vi,,) 



(62) 



Out 



Vin + Jn 



Ti is solution of Jmax,tTi + 
Tl = T3 

T2 = T4 = T5 = Tg = T7 = 



2V/„Ti -L = 



(63) 
(64) 
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Figure 12: Cases 1 and 2 
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Case 3: Vout is reached by reaching AmaxA- Vp is not reached. Vp is reached by reaching 
Ti = T3 = Ts = r7 = Amax,i/ Jmax,i Vp, Ti and tq are solutions of: 



Vout 

L 



A. 



1 



+ V'-A 



2Jrr 



max.i ' ( 



(65) 



Case 4: Vout is reached without reaching Amax,i- Vp is not reached. Vp is reached by reaching 
Amax,i- Ti = T3 = Amax,i/ Jmax.i and Vp, T2, T5 are solutions of: 

Al 



vi. 

Vout 

L 



^ + Vin + A 

V' -t'^ J 
F 5 '^rnax.i 

' 5 •^max,z I 



+ 



+T2T^/„ + 2r5T/; 



(66) 



Case 5: Vout is reached by reaching Amax.i- Vp is not reached. Vp is reached without reaching 
Amax,i- T5 = T7 = ^max.i/^max.i, T2 = 0, Ti = T3, Ti , and Vp arc solutions of: 



Vout 

L 



Vir, 



J, 



max.! ' 1 



iT6 



2J^ 



I ^ ' 1 '^max,i 

+ (4tiV^/„ + 2t6V^; - ^ma.,»T|) 

-2A^ ) 

max,i J 



(67) 
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Case 6: Vout is reached without reaching Amax,i- Vp is not reached. Vp is reached without reaching 
Amax,i- Ti = T3, Ts = Ty, Ti , T5 and are solutions of (figure 13) : 

Vp — T^JmaxA + 2riV7„ 

L = Jmax,i + TfJmax,i + 2tiV7„ + 2TfVp 



Figure 13: Case 6 
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• Case 7: Vout is reached by reaching Amax.i- Vp is reached. Vp is reached by reaching A„ 



\VF~Vln\ > 
\Vout-Vp\ > 

L > 



Jm.a 



{2J„iax,iVp + 2A^„^ j) Vin 
'Jmax.iVp + A^^^^jVp 



(69) 



max, I ' 

2v4? 



{^Jmax^i^Oi 
Jmax.i ^ Out 



Vp 



Tl = T3 = T5 = = ^ 



T2 



(70) 



• Case 8: Vout is reached without reaching Amax,i- Vp is reached. Vp is reached by reaching Ajnax,i- 



^ . (2./™„x,,Vf+2A^„^ ;)Vj„ + J„„x,,V^+A^ 
^ ^ T 



(71) 



V J„.ax Vo^f (4 V> + Vo^t )+A„„x 
2J™„x 



n = T3 = j^^^-^ 



T5 = T7 : 

I re = 



Jrr 



Vp-Vout 



(72) 
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• Case 9: Vout is reached by reaching Amax,i- Vf is reached. Vp is reached without reaching A„ 

\Vf-ViJ < 



\Vout-VF\ > 



(73) 



n = T3 = 

r2 = 
T5 = T7 = 



(74) 



• Case 10: Vout is reached without reaching Amax,i- Vp is reached. Vp is reached without reaching 
Amax.i (figure 14) : 



|1^F-V/„| < 

\Vout - Vf\ < 4ff7i 



L > 



VJma^VF{4:Vl„ + VF)+A„,^:,Vl„ 
VJma^V0uti4:VF+V0^t)+A^a^VF 



(75) 



n = Ta 
r2 =0 

T5 

I T6 =0 





-Vf 















(76) 



Figure 14: Case 10 
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In the simulations, equation 2 is resolved by the Cardan method. The systems of equation 65, 66, 67 
and 68 are non-linear. They are resolved numerically using the Newton-Raphson method. 

To use this algorithm, take a segment of 0.01m to be covered. Take Vf„ — Q.2m/s, Vout — O.lm/s, 
Vp — 0.5m/ s. The machine parameters are those in table 1. 

A^ / 

\Vout — Vin\ < j"""''' ; thus there will be no phase 2. L> {Vout + Vin) J ^^y^^'^ ■ the length of the 
segment thus allows Vout to be reached. 
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\VF-Vin\ < ^ffff and \Vout-VF\ < ^fff and L < VJrr.a.Vo.t{4V^+Vo.t)+A^a.v^ r^^^ ^^^^^^ ^^^^ 
therefore not allow Vp to be reached and the differences in feed rates will not allow maximum acceleration 
to be reached. This means case 6 applies. This gives = 0.4418m/s, 

n = T3 = 0.07776s. T5=Tr = 0.09245s. 
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